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We present a new calibration of the elemental-abundance data for Asteroid 433 Eros taken by the X-ray 
spectrometer (XRS) aboard the NEAR-Shoemaker spacecraft. (NEAR is an acronym for "Near-Earth Asteroid 
Rendezvous.") Quantification of the asteroid surface elemental abundance ratios depends critically on 
accurate knowledge of the incident solar X-ray spectrum, which was monitored simuitaneously with 
asteroid observations. Previously published results suffered from incompletely characterized systematic 
uncertainties due to an imperfect ground calibration of the NEAR gas solar monitor. The solar monitor 
response function and associated uncertainties have now been characterized by cross-calibration of a 
large sample of NEAR soiar monitor flight data against contemporary broadband solar X-ray data from the 
Earth-orbiting GOES-8 (Geostationary Operational Environmental Satellite). The results have been used to 
analyze XRS spectra acquired from Eros during eight major soiar flares (including three that have not 
previously been reported). The end product of this analysis is a revised set of Eros surface elemental 
abundance ratios with new error estimates that more accurately reflect the remaining uncertainties in 
the soiar flare spectra: Mg/Si = 0.753 +0.07S/-0. 055, Al/Si = 0.069 ±0.055, S/Si =0.005 ± 0.008. Ca/Si = 
0,060 + 0.023/— 0.024, and Fe/Si = 1.678 + 0.338/— 0.320. These revised abundance ratios are consistent 
within cited uncertainties with the results of Nittler et al. (Nittler, L,R„ and 14 colleagues, 2001. Meteorit. 
Planet. Sci. 36, 1673-1695] and thus support the prior conclusions that 433 Eros has a major-element 
composition similar to ordinary chondrites with the exception of a strong depletion in sulfur, most likely 
caused by space weathering. 

Published by Elsevier Inc. 


1. Introduction 

The Near Earth Asteroid Rendezvous (NEAR)-Shoemaker space- 
craft orbited the Asteroid 433 Eros between February 2000 and 
February 2001. Among its mission objectives {Cheng, 1997) was 
the characterization of the surface composition of Eros using X-ray 
(XRS, Goldsten et ai„ 1997) and gamma-ray spectrometers (GRS, 
Goidsten et al, 1997). XRS results were published during and im- 
mediately after Che mission by Trombka et ai. (2000) and .Nittler et 
at. (2001), but were plagued by systematic uncertainties stemming 
from the unlucky failure of one solar monitor and the incomplete 
pre-launch calibration of the other. This paper describes the sub- 
sequent effort at post-launch calibration of the NEAR gas solar 
monitor and the resulting improvement in our understanding of 
the XRS-derived elemental composition of 433 Eros. 

The X-ray fluorescence emission measured by NEAR is a prod- 
uct of the interaction of soiar X-rays with the asteroid's surface. As 
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the asteroid is exposed to X-rays from the upper soiar atmosphere, 
a fraction of the incident X-rays is photoelectrically absorbed by 
various inner electron shells of atoms in the surface, resulting in 
the emission of either an Auger electron or a fluorescent X-ray, 
with energies characteristic (Moseley, 1913) of the atoms that pro- 
duce them. 

Each electron shell can only be induced to fluoresce by that 
portion of the soiar X-ray spectrum that fails above its binding 
energy. Moreover, fluorescent X-rays emitted by one atom may 
be absorbed by another atom in the asteroid's surface prior to 
escape into space. The resulting spectrum of fluorescent lines de- 
pends both on the abundance of each element in the asteroid’s 
surface and on the spectral distribution of the incident X-rays. Ac- 
curate knowledge of the incident solar spectrum is thus critical to 
converting the spectra received by the NEAR XRS into elemental 
abundances on the asteroid's surface. 

NEAR, therefore, carried two X-ray spectrometers dedicated to 
monitoring the solar spectrum between t and 10 keV, One of 
these, a Si PIN photodiode, failed prior to orbit insertion. The sec- 
ond consisted of a gas-fiited proportional counter similar to the 
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Fig, 1, Theoretical (CHIANTI 5.2; Dere et aL 1997; Landi er ai.. 2006) single-remperature solar spectra at three plasma temperatures. Vertical dashed lines represent the 
K edge energies of major elements measured by the NEAR XR5. 


asteroid-pointing detectors, but covered by a graded filter (Clark 
et a!,. 1995} designed to attenuate the incoming soiar flux and 
increase the dynamic range of the detector. Unfortunately, the re- 
sponse of this filter, particularly at off-normal solar incidence an- 
gles, was not adequately determined in the laboratory prior to 
launch. 

This paper describes the post-launch calibration of the NEAR 
gas soiar monitor and its application to the analysis of XRS spectra 
of 433 Eros acquired during eight major soiar flares. The response 
of the graded filter and the uncertainties therein were modeled 
based on measurements and tolerances as recorded in the engi- 
neering diagrams. Resulting models of the gas solar monitor’s re- 
sponse function were tested against flight data and checked against 
contemporaneous broadband X-ray data from the Earth-orbiting 
GOES-8 (Geostationary Operational Environmental Satellite) to re- 
duce the range of possible solar monitor response functions to 
those that were consistent with the available data. Finally, the 
improved soiar results were used to analyze data from eight ma- 
jor solar flares from which asteroid fluorescence was recorded 
by NEAR, resulting, in a recalibrated elemental composition for 
433 Eros. The remaining uncertainties in the solar spectra were 
incorporated into the error analysis for the asteroid’s composition. 

2. The soiar X~ray spectrum 

The solar X-ray output has been observed by instruments 
aboard a number of Earth-orbiting satellites. These have included 
various GOES (Geostationary Operation Environmental Satellites) 
broadband detectors; the SMM Bent Crystal Spectrometer (Acton et 
al„ 1980) in the 1980s (1980, then 1984-1989); and the YOHKOH 
BCS and SXT (Bragg Crystal Spectrometer and Soft X-ray Telescope; 
Yoshimori et al„ 1991 } in the 1990s. After the NEAR mission, solar 
X-ray observations have also been conducted by the CORGNAS-F 
X-ray spectrometer (RES1K; Sylwester et al„ 2005) operated from 
2001-2003, and RHESSi (Reuven Ramaty High-Energy Solar Spec- 
troscopic Imager; Lin er al, 2002), launched in 2002. 

The solar X-ray spectrum consists of a set of emission lines su- 
perimposed on a continuum. Emission lines in the region observed 
by NEAR (1-10 keV) come from highly ionized aroms (H and He- 


like electron shells). The continuum is produced by bremsstrahlung 
emission, primarily from hydrogen, and declines steeply with in- 
creasing energy. Both the slope of the bremsstrahlung continuum 
and the intensities of the various emission lines vary greatly with 
electron temperature (Fig. 1). Electron temperatures (T e ) in the 
corona generally are in the range --2-3G x 10 6 K (e.g., Phillips, 
2004), with the highest temperatures being reached only during 
major solar flares. 

The amount of material at a given T e is described by the “emis- 
sion measure" (EM): 

EMs jn 2 dV, (1) 

where n represents the plasma density and V is the volume of the 
hot plasma. The total spectrum for an isothermal plasma, then, is: 

Remitted f'e) = (Eff(2., TV } 4- £fb(4, 7 e ) } j uj. dV (2 j 

(Garcia, 1994), where Eq, is the emission line spectrum and Eq is 
the bremsstrahlung continuum. 

Generally, however, the solar spectrum is multithermal (e.g., 
Feldman et al„ 1995; McTieman et al„ 1999). The thermal portion 
of the solar X-ray spectrum, then, can be represented by the inte- 
gral over ail temperatures of the amount of plasma at a given tem- 
perature (known as the "differential emission measure," DEM(T)) 
multiplied by the spectrum appropriate to that temperature: 

&ctai(X) = j T) x DEM (J}dT. (3) 

Typically, DEM(T) is calculated by observing the relative intensities 
of lines in the solar spectrum, since different species of ions have 
different emitting efficiencies as functions of temperature. 

Because the NEAR gas soiar monitor could not resolve individ- 
ual spectral lines, its spectra do not contain enough information 
to permit full inversions of DEM(T). The typical structure of the 
DEMfT) of a solar flare represents two plasma components, one 
emitting at 5-10 MK and a hotter component emitting at 16- 
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Fig, 2 . Singie-temperature (top) vs twa-temperature fit to a NEAR solar monitor 
spectrum from 17-Jtme-2GGG. The smooth dashed traces are the models, while the 
noisy solid ones are the NEAR data. In the lawef plot, the dash-dotted traces are 
the two isothermal components of the fit. Two-temperature modeling was used 
throughout the present analysis. 


25 MK (e.g., Antonucci and Dodero, 1995; Kepa et al„ 2006'.' We 
have therefore made the simplifying assumption that the solar X- 
ray spectrum can be adequately described by a two-temperature 
mode!, so that; 

<!>»,oxi\Q-) = EMi x^(A,r,) + EM 2 x^(>„r2). (4) 

Fig. 2 illustrates single-temperature and two-temperature fits to 
a NEAR gas solar monitor flare spectrum. 

2:1. Elemental abundances in the solar atmosphere 

The line emission in the solar X-ray spectrum depends on the 
abundances of the emitting elements in the solar atmosphere as 
well as on the electron temperature. Generally, the X-ray emit- 
ting plasmas of solar flares reflect the composition of the solar 
corona, which differs substantially from that of the photosphere; 
in the corona, elemenrs (other than hydrogen) with high first ion- 
ization potentials are depleted relative to elements with low first 
ionization potentials (the ‘'FIP” effect). Although this relative frac- 
tionation according to FiP is well established, the absolute nor- 
malization of corona] abundances relative to photospheric ones is 
somewhat controversial (e.g., Fiudra and Schmeiz, 1999). !n this 
work, we consider the published coronal abundances of Meyer 
(1985), based on observations in the extreme ultraviolet and X-ray 
ranges, and of Feldman et al, (1992). 

Synthetic solar spectra were generated using the coronal abun- 
dances of both Meyer (1985) and Feldman et al, (1992) as inputs to 


f i he solar flares examined in these papers represented GOES classes C2-X2 aril 
C8.1, respectively. 
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the CHIANTI 5,2 code (Dere et al„ 1997; Land! et al„ 2006). In the 
latter calculation, a number of minor elements not included in the 
work of Feldman et al. (1992) were assigned abundances based on 
the photospheric abundances of Grevesse and Sauvai (1998), with 
elements with low first-ionization potentials (FiP) increased by a 
factor of 3.5 for the "FIP" effect. These are the same values used 
in the “sun^coronaLexf file included with the CHIANT! Soiarsoft 
package (Land! et al., 2006). 

Neither abundance set was entirely successful in modeling the 
NEAR gas solar monitor data: the spectra generated with the 
Feldman et ai. abundances produced greatly superior matches in 
the region below 2.5 keV, but performed poorly above 6 keV; 
whereas the models using the Meyer abundances produced excel- 
lent fits to the data above 6 keV. Since the hydrogen-dominated 
bremsstrahlung continuum is an important contributor to the 
X-ray spectrum in the high-energy region where emission lines 
are relatively sparse, the overall metailicity of the Feldman et a!. 
(1992) model was adjusted downward so that the Sr/H ratio would 
be the same as that of Meyer (1985): Si/H = Fe/H = 3.9 x 10 -s . 
The resulting mode! produced a good fit to the NEAR data at all 
energies. Subsequently, additional models were tested with coro- 
nal metallidties varied to be slightly higher and slightly lower than 
those of the Meyer (1985) model. The two most successful mod- 
els were those with Si/H = 3.5 x 10~ 5 (“Coronal 1“ in Table 1) and 
Si/H = 3.9 x IG -3 (“Coronal 2"). Abundances in all relevant models 
are included in Table 1. 

Fiare-to-flare variation in solar elemental abundances has been 
observed by several investigators (Sylwester et al., 1984; Sterling 
et ai., 1993; Fiudra and Schmeiz, 1999), For example, Sylwester 
et ai. (1998) found that the calcium abundance varies from flare 
to flare over a factor of about 3.5 and is uncorrelated with elec- 
tron temperature or GOES classification. Within a given flare and 
among flares associated with the same active region, however, the 
abundance was found to change only slowly. The NEAR gas so- 
lar monitor data from the major solar flare of 27 December 2000 
(Section 6.5) proved impossible to fit with the same abundance 
sets that had worked for the other NEAR flare data. In particu- 
lar, the 6.5 keV iron structure was clearly weaker with respect to 
the surrounding continuum than was the case in the spectra of 
other solar flares. For this reason, additional models were applied 
to this flare in which the abundance of iron was reduced by vary- 
ing amounts. 

3. The NEAR gas proportional counters 

The energy resolution of the NEAR gas proportional counters, 
including the gas solar monitor and the three asteroid-pointing de- 
tectors, varied with energy according to FWHM = 0.340£ 1/2 (Starr 
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Fig, 3. Geometry of the NEAR gas solar monitor graded filter, left: Engineering diagram of the graded filter assembly as mounted above the NEAR gas solar monitor tube. 
The filter pieces {items 4-T) each measured 2 x 3.5 inches. Right: Schematic cross-section of the graded filter pinhole with nominal dimensions. The gas solar monitor was 
mounted on the sunward-facing deck of the NEAR-Shoemaker spacecraft so that the sensor face was parallel with the solar panels. Thus, the angle marked is the same as 
the “panel angie :1 of the spacecraft. 


et al„ 2000), where FWHM is the fuil width at half maximum of 
a Gaussian X-ray line. As discussed earlier, this resolution is not 
sufficient to resolve individual spectra lines in the solar X-ray spec- 
trum. precluding inversion of the full DEM(T) (Eq. (3)) and limiting 
the information that can be obtained on solar elemental abun- 
dances. A potential drawback to our method of fitting solar moni- 
tor spectra to a two-temperature model (Fig. 2) is that the further 
the actual solar DEM{T) departs from the sum of two delta func- 
tions at two temperatures {for example, if there is a major mass 
of plasma contributing at a third temperature, or if the plasma 
temperature varies over a wide range such that a roughly equal 
amount of X-ray flux is contributed from each temperature incre- 
ment AT), the poorer a representation of the actual solar spectrum 
our model will be. 

Additional complications to modeling spectra from the NEAR 
gas detectors are due to escape peaks and pulse pile-up effects. 
The NEAR proportions) counters were filled with P-10 {Starr et 
ai„ 2000), a standard mixture of 90% Ar and 10% G-L. Methane 
does not significantly absorb X-rays in the 1-10 keV range; its 
purpose is to absorb the low-energy photons emitted by the ar- 
gon in order to prevent spurious pulses. Argon, however, can cause 
an escape peak: if an incident X-ray induces fluorescence from an 
Ar atom and the fluorescent X-ray escapes the detector, the en- 
ergy measured will be 2.97 keV (the K a energy for Ar) below the 
original energy. Thus, an "escape peak" appears in the output spec- 
trum for all incident energies higher than the binding energy of 
Ar (3.2 keV). For the NEAR detectors, the escape fraction was cal- 


culated using Monte Carlo methods (R. Starr, 8 December 2004, 
personal communication), and was found to be well described by 
a linear function of energy. 

At high count rates, such as those that occurred in the gas so- 
lar monitor during major solar flares, pulse addition in the detector 
electronics caused a pile-up effect in the entire spectrum, such that 
each channel would represent a higher energy than it did at low 
count rates. The size of this “zero shift" in the energy vs chan- 
nel calibration increased with count rate, and was treated as a free 
parameter in the spectral fitting. For context, during the flares de- 
scribed in the present work, the overall calibration shifts ranged 
from ss2-56 eV. 

4. Modeling the NEAR gas solar monitor 

The NEAR gas solar monitor was covered by a graded filter 
(Clark et ai., 1995) intended to attenuate the incoming solar flux 
and increase the dynamic range of the detector (Fig. 3). Since the 
Sun’s X-ray emission decreases steeply with increasing energy, the 
filter was designed to admit a greater proportion of “high-energy” 
(=«4-9 keV) X-rays than low-energy ones. 

The graded filter was constructed from thin layers of beryl- 
lium. aluminum, and the plastics Delrin and Kapton (Fig. 3). To 
prevent the proportional counter from being saturated, most of its 
area was covered by the Al layer, A broad slot (0.295 ± 0.005 by 
1.181 ± 0.005 indies, or 2,25 cm^) was cut out of the Ai layer 
to admit a portion of the original incident solar radiation. The 
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Table 2 

Nominal dimensions and engineering tolerances for the graded filter. 
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’ As measured by Goidsten (2003, personal communication). Tolerance not spec- 
ified in engineering drawings. 

h Not regarded as a critical dimension according to the engineering drawings: 
uncertainties unknown. 

c According to Goidsten (2003), this "is the 100 HN process of Kapton by Dupont. 
The single point thickness is 0.S5-1.15 mils for a 1.0 mil nominal thickness." 

Window of the off-the-shelf gas counter. Not part of the graded filter. 

Delrin and Be layers underneath were pierced by concentric pin- 
holes intended to allow a small Fraction of the incident low-energy 
X-rays to reach the proportional counter underneath. Meanwhile, 
higher-energy X-rays would be relatively unimpeded by these ma- 
terials. 

The nominal dimensions of the pinholes and thicknesses of the 
graded filter layers are summarized in Table 2. These dimensions 
are critical for accurately determining the response of the gas solar 
monitor. Unfortunately, variations in these dimensions, even within 
the engineering tolerances of the filter, have substantial effects on 
the energy and angular response of the detector. The specified 
measurements come from the engineering drawings tor the graded 
filter (Goidsten, 2003, personal communication). Given pinhole di- 
ameter uncertainties are the nominal engineering tolerances and 
material thicknesses are the specified stock thicknesses from the 
manufacturers. 

Because the Delrin pinhole diameter was similar in dimension 
to the Delrin stock thickness (Pig. 3), the stock thickness is crit- 
ically influential in determining the response function of filtered 
solar monitor at off-normal solar incidence angles. Unlike the pin- 
hole diameters, however, the thicknesses of the filter layers were 
not marked as critical dimensions in the engineering diagrams. 
Regrettably, it has proven impossible to recover any pre-iaunch 
measurements that may have been made of these thicknesses. 
However, a leftover sample of the Be filter foil was later recov- 
ered by Goidsten et at. ( 1997) and measured after the fact to be on 
average 0.0042 incites thick, with individual measurements rang- 
ing from 0.0039 to 0.0046 inches. The Kapton stock was stated 
by the manufacturer (Richmond Aircraft Products) to be t.O mil 
(0.001 in.) in thickness, with single-point thicknesses ranging from 
0.85 to 1.15 mils. No estimates of uncertainty were available on the 
thicknesses of the Delrin or aluminum stock or on the densities of 
the Delrin and Kapton plastics. 

Attenuation coefficients for materials used in the proportional 
counter and graded filter were obtained from NIST (Hubbeil, 1997, 
available at http://www.physics.nist.gov/PhysRefData/XrayMassCoef 
/cover.html). 

4.1. Post-launch calibration of the filtered gas solar monitor 

In order to reduce the uncertainties associated with the limi- 
tations of the pre-iaunch graded filter calibration, flight data from 
the NEAR mission were used to determine the range of plausible 
graded filter parameters consistent with the solar spectra collected 
during the mission. The parameter space within the engineering 
tolerances of the graded filter was searched systematically. Each 


filter mode! was tested against actual gas solar monitor output 
at a wide range of incident angles and solar activity levels. One 
hundred sixty NEAR spectra were selected from solar flares that 
occurred between 4 May 2000 and 1 July 2000. Most of these 
flares occurred when the sunlit side of the asteroid was not in 
the field of view of the NEAR XR5, so are not represented in tire 
asteroid-pointing data. For each NEAR spectrum, which represents 
a 50-s integration of the solar monitor, we generated a correspond- 
ing set of synthetic gas solar monitor spectra. Theoretical (CHIANTI 
5.2: Dere et a!., 1997; Landi et a!., 2006) models of solar output 
at thirty-nine plasma temperatures between 2 and 40 MK were 
put through the response function of the candidate graded fil- 
ter model, calculated at the correct solar incidence angle for the 
time of the observed spectrum, and then convolved with the re- 
sponse function of the gas counter. The best two-temperature fit 
according to a x 2 test was then recorded, along with the tem- 
peratures and emission measures of the two isothermal compo- 
nents. 

The free parameters in the fitting process were the tempera- 
tures of the two single-temperature spectra to be combined; the 
relative proportions of the law- and high-temperature compo- 
nents; the overall normalization, which is related to the emission 
measure of the spectrum; and the zero shift of the detector due to 
pile-up, as described in Section 3. Over a targe sample of incidence 
angles and solar spectra, the lowest average x 2 value should cor- 
respond to the graded-fiiter model that best described the flight 
hardware. 

In addition to the goodness of the fit to the solar monitor out- 
put itself, the second test of our filter models was provided by 
contemporary data on solar X-ray emission provided by the GOES- 
8 satellite. Throughout June and August 2000, the asteroid-Sun- 
Earth angle was below 30°. so that during this period (but not, for 
example, during December 2000) GOES and NEAR were observing 
substantially the same portion of the Sun. 

GOES-8 was equipped with two broadband X-ray channels: a 
short-wavelength channel with a bandpass of 0.5 to 4 A (24.8 
to 3.1 keV) and a long-wavelength channel at 1 to 8 A (12.4 to 
1.5 keV), Several investigators (Garcia, 1994; Feldman et a!„ 1995; 
Thomas et al„ 1985) have published formulae for deriving tem- 
peratures from one or both of the GOES channel outputs. These 
were used by Trombka et ai. (2000) in the earliest NEAR XRS 
paper to interpret the asteroid-pointing data, in place of the so- 
lar monitor; and similarly by Nittler et a!. (2001 ) to interpret the 
“quiet sun” (non-flare) data. However, all of these GOES models 
are based on an isothermal approximation of the flare plasma, 
in order to test our new muititherma! fits to the solar moni- 
tor data, we instead put each synthetic spectrum directly through 
the two GOES transfer functions (Thomas, 2003; Schwartz, 1998; 
Garcia, 1994) to model the outputs of the two channels. The ratio 
of the modeled channel outputs was then compared to the ratio 
of tile actual GOES channel outputs (“GOES ratio”) from the same 
time period, corrected for light time. 

Based on the results of the x 2 fitting and comparison to GOES 
data, eleven graded-fiiter models (Table 3) were selected for use 
in the asteroid-fluorescence data analysis. Eight of these represent 
the best y 2 matches to the solar monitor data on the entire 160- 
spectrum test sample and the best matches to the GOES channel 
ratios in the subset (56 of the 160) taken from June 21 through 
July 1, when the Earth-Sun- asteroid angle was <16.5“. With these 
eight models, the median y 1 values for the 160 fits in the test 
sample were under 1.765 and median RMS matches of the calcu- 
lated vs the measured GOES ratios under 0.032, 

Three additional models were included. Model 9 was selected 
because it produced the lowest median x 2 fit to the solar monitor 
data of any that was tested (yfi ~ 1.658), although its performance 
on the GOES ratio test was slightly worse than that of a number of 
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Table 3 

Selected graded filter models 
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Table 4 

Geometric information for eight solar flares. 



other models. (Model 1, which produced the lowest median RMS 
to the GOES channel ratios— 0.022— was already included among 
the first eight models.) Models 10 and 11 represented "local min- 
ima" in both x 2 and GOES ratio RMS comparisons; 10 was the 
only model that performed well in which the Kapton layer was 
as thin or thinner than its nominal stock thickness, and 11 was 
the best-performing model among those in which the Deirin was 
as thick or thicker than its nominal stock thickness. Although the 
eight overall best models ail had thick Kapton and thin Deirin 
layers, the performance of these two models was dose enough 
that they could not be ruled out on the basis of the available 
data. 


5. Application of the solar monitor results to the determination 
of asteroid elemental abundance ratios 

We have applied the new solar monitor models in a new analy- 
sis of asteroid XRS spectra for eight solar flares; five from May-July 
2000, two from December 2000, and one from 2 January 2001. 
Note that these include the five flares analyzed by Nittler et al. 
(2001) and three additional flares for which data have not previ- 
ously been reported. The time and geometric information for the 
analyzed flares are provided in Table 4. 

The procedure used to determine surface abundance ratios of 
Eros from NEAR XRS data acquired during solar flares followed 
closely that used by Nittler et at. (2001). The steps to this pro- 
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Fig> 4, Assembly of the NEAR X-ray asteroid-pointing sensors (Fig. 7 from Gofdsien 
et al. (1997), reprinted by kind permission of Springer Science and Business Media). 
To enabie analysis of the low-energy X-ray lines, one detector was filtered with a 
thin layer of A] and another was similarly equipped with a Mg filrer. The center 
detector remained unfiltered. The rectangular faces of the gas counter tubes (“X-ray 
detectors 1 ') were each 2 x 3.5 inches, as with the gas solar monitor tube (Fig. 3). 
The same model tubes were used in both cases. 


cess were: determination of the appropriate solar X-ray spectrum 
from solar monitor data, extraction of the photon flux in each of 
the major fluorescent lines from the asteroid-pointing spectra, gen- 
eration of calibration curves relating photon ratios to abundance 
ratios for a specific incident solar spectrum and viewing geometry, 
and conversion of photon flux ratios to abundance ratios. 

Solar models for each flare were produced using all eleven 
graded filter models, selected as described above. For each fil- 
ter model, each 50-s integration of the gas solar monitor was fit 
with a two-temperature model of the solar spectrum, resulting in 
a model solar spectrum with sufficient energy resolution to dis- 
tinguish individual solar emission lines. The 50-s two-component 
mode! spectra were then co-added to produce a time-integrated 
solar flare spectrum representing the same time interval as the co- 
added asteroid spectra (Section 5,2). 

The fitting process used to extract asteroid fluorescence counts 
from the asteroid-pointing detector spectra is identical to that de- 
scribed in Nittier et al. (2001 ). with the exception that the on- 
board energy calibration data has now been analyzed and applied 
to constrain the properties of the proportional counters, as de- 
scribed below. 

5,1. Calibration of the asteroid-pointing detectors with the on-board 
^resources 

Nittier et al. (2001) allowed the gains and zeroes of the 
three asteroid-pointing detectors to float as free parameters in 
their previous analysis of data from five NEAR solar flares. How- 
ever, the NEAR XRS included three radioactive 55 Fe calibration 
sources (“Calibration Rod Assembly" in Fig. 4) to monitor the gain 
(eV/channel), zero (energy of channel zero), and energy resolu- 
tion of the asteroid-pointing detectors. Several times throughout 


Table 5 

XRS asteroid-pointing detector calibration data fitting results. Detector resolution is 


a function of energy: FWKM(E) — FWHM parameter x 
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Table 6 

Unfiltered detector characteristics derived from flare data. 
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the mission (Table 5) the calibration sources were rotated into the 
fields of view of their respective detectors for the acquisition of 
calibration data. Making use of the in-flight calibration data re- 
duces the number of free parameters, increasing confidence in the 
fitting results. This is especially important for the two filtered de- 
tectors, in which count rates from the asteroid were always low in 
comparison with the cosmic-ray-iuduced background. 

The S5 Fe sources emitted X-rays in the Mn (88%) and K/j 
(12%) lines at 5.898 and 6.490 keV, respectively. The detectors also 
registered escape peaks for the two lines at 2.941 and 3,533 keV. 
For every available set of calibration data, the detectors’ response 
to the pair of monoenergetic lines received from the 55 Fe source 
was modeled and fit to the calibration data. Gaussian shapes were 
assumed for all four peaks. The scale of the cosmic-ray back- 
ground, whose shape was established by calibration-source- and 
asteroid-free data taken within 24 h of the calibration data, was 
also allowed to float, but usually remained within 1% (and always 
within 2%, for the two filtered detectors) of simple scaling by in- 
tegration time. Calibration source fitting results are summarized in 
Table 5. 

The 55 Fe source for the unfiitered detector was much weaker 
than the sources for the other two detectors, resulting in poor 
statistics on its calibration data and consequently relatively large 
error bars on its gain, zero, and FWHM parameters in Table 5. 
However, the unfiltered detector had the best statistics of the three 
in the asteroid fluorescence during the solar flares themselves, 
and more precise detector characteristics could be derived from 
the flare spectra themselves (Table 6). In contrast, the calibration 
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Fig. 5. XRS da ta from the three asteroid-pointing detectors. Left to right: Raw data from the solar Hare of 27-Dec-2G0G; background-subtracted data from 27-l}ec-20QO; 
background-subtracted data from the flare of 17-Jun-2000; background-subtracted data from the flare of 28-Dec- 2000. The step at channel 75 in the raw data is the rise-time 
discriminator limit. Solid lines represent the best Gaussian fits to individual XRF peaks; dashed traces represent the sum of ail components. From left to right, the peaks are 
Mg, A!, Si, S, Ca, and Fe {K,, and K,)). 


source data gave better statistics than the flare data for the two 
filtered detectors. 

The energy resolution was the most variable of the detector 
characteristics, becoming degraded whenever the instrument was 
shut off, for example when the spacecraft was placed in safe 
mode, and taking several weeks to approach its best value (FVVHM 
parameter = 0.340). This can be observed in Table 5, as the resolu- 
tions of the AI and Mg detectors improve slowly from April to late 
May. 

5.2. Extracting photon counts from the asteroid-pointing detector dam 

As with the solar monitor data, the spectra from the asteroid- 
pointing detectors were recorded in 50-$ integrations. For each so- 
lar flare, all the spectral integrations from that flare were co-added 
to generate a single summed flare spectrum for each asteroid- 
pointing detector. Cosmic-ray induced background, present in ail 
of the gas proportional counters, was then subtracted from each 
co-added spectrum prior to fitting. Since this background was con- 
stant on time scales of a few hours, background spectra were de- 
rived for each detector from data taken within 24 h of the flare 
while the asteroid was not in the XRS field of view. These "dark" 
spectra were added together and scaled to the live time of the rel- 
evant flare spectra. Examples of co-added asteroid spectra before 
and after background subtraction are plotted for the 21 December 
2000 flare in Fig. 5. 

Following background removal, the spectra from the three de- 
tectors were fit simultaneously, allowing the Mg- and Al-filters to 
perform their design function: to compensate for the low resolu- 
tion of the proportional counters at the K a energies of the light 
elements Mg, Ai, and Si, The free parameters in the spectral fit- 
ting were the intensities of the fluorescence lines from each major 
element (Mg, Ai, Si, S, Ca, and Fe) entering the top of the instru- 
ment, as well as the gains, zeroes and FVVHM parameters of the 
detectors. For the unfiltered detector, the detector parameters were 
allowed to vary freely: while for the Mg- and Al-filtered detec- 


tors, they were allowed to float within the limits derived from 
the calibration-source data (Table 5) closest to the date of the 
flare. 

The fitting process for three solar flares is illustrated in Fig. 5. 
In each plot, the grey lines with error bars represent the data from 
each asteroid-pointing detector. Error bars are la, reflecting count- 
ing statistics. The leftmost plot shows the raw asteroid-pointing 
data from the flare of 27 December 2000, whereas the three re- 
maining plots show the data after the cosmic-ray background has 
been subtracted. The fit to each of the individual fluorescence lines 
is represented by a solid Caussian trace. From lowest- to highest- 
energy, these are the K a lines of Mg, AI, Si, S, Ca, and Fe and the 
K ft line of Fe. The effect of the balanced -filter system can easily 
be observed in these plots: the Si line, which is the most promi- 
nent feature in the unfiltered data (top plots), has been suppressed 
in the data from the two filtered detectors. The superimposed 
dashed trace is the overall fit to the data. Similar plots represent- 
ing the five previously analyzed flares were printed in Niftier et ai. 
( 2001 ). 

5.3. Conversion of photon ratios to asteroid elemental ratios 

The amount of fluorescence in any line that escapes the asteroid 
depends not only on the abundance of the fluorescing element, but 
also on the abundance of all other elements in the asteroid surface 
that can absorb incident solar radiation on the way in or fluores- 
cent radiation on the way out. Because of this, it is necessary to 
make some assumptions about the entire composition of the aster- 
oid in order to interpret the XRS data, including elements whose 
fluorescence could not be detected. We based our determination of 
the relationships between fluorescence ratios and elemental ratios 
on data from meteorites. 

Once the incident solar spectrum was modeled, theoretical 
X-ray spectra (including both fluorescence and coherent scattering 
components) that would be induced by the model solar spectrum 
at the applicable incidence, emission, and phase angles were cal- 
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Table 7 

Meteorite samples used to convert photon ratios to element ratios. 
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cutated for a sample of meteorite compositions of widely different 
classes (Table 7), This is the same method described and used 
by Nittler et al. (2001). The model X-ray spectra for the differ- 
ent meteorite compositions were then fit using exactly the same 
method as had been applied to the actual XRS spectra to derive 
model fluorescence ratios. Calibration curves were calculated by 
plotting actual meteorite elemental abundance ratios against mod- 
eled meteorite fluorescence ratios and fitting quadratic curves to 
these points. Each quadratic fit provided the formula for convert- 
ing fluorescent- photon ratios to elemental ratios for a particular 
incident solar spectrum and geometry. Fig. 6 illustrates the cali- 
bration curves for the solar spectrum and geometry of the major 
flare from 17 June 2000. 

As in Nittler et al. (2001 ), the coherent scattering background 
has been accounted for by inclusion in these calculated meteorite 
fluorescence ratios rather than by any attempt to subtract it from 
the fluorescent photon signal. The induced A! fluorescence from 
the Al filter in front of the Al-filtered detector has been treated 
similarly. Because the theoretical spectra are treated exactly the 
same way as the real data, this does not introduce significant 
errors. However, because the derived fluorescence counts include 
scattered photons and filter fluorescence, the photon ratios for 
Al/Si and S/Si in the calibration curves (Fig. 6) are nonzero even 
at zero concentrations of A! and S. 

The angles used in these calculations were averages over those 
provided in the XRS data products for the individual spectral in- 
tegrations during a flare. The phase angles were exact, but the 
incidence and emission angles were averages over all the "plates" 
of the shape mode! derived from the NEAR laser rangefinder data 
(Zuber et a!., 2000} that were illuminated and within the XRS 
field of view at the time. These derived angles varied much more 
during flares than did the phase angles (Table 4) and thus, in prin- 
ciple, calibration curves should be determined for each integration. 
However, the statistics did not permit separate fitting of integra- 
tions over small time portions within a single flare. Modeling of 
expected X-ray spectra for an ordinary chondrite composition illu- 


minated by a typical solar flare at a fixed phase angle, but with 
varying incidence and emission angles, indicated that the uncer- 
tainty introduced by using averaged incidence and emission angles 
is significantly smaller than that caused by other sources of er- 
ror in our analysis, for example the solar modeling. Moreover, the 
incidence and emission angles that matter physically for X-ray in- 
teractions are on a much smaller scale than any possible shape 
model of the asteroid, so this issue is a fundamental limitation in 
the analysis of NEAR XRS data. 

In their XRS analysis, Nittler et al. (2001) applied an additional 
correction to the data to take into account the fact that in mete- 
orites (and hence probably stony asteroids), the elements are not 
spread uniformly throughout the asteroidal material, but rather 
tend to be segregated into different minerals on spatial scales 
that are important to X-ray absorption and scattering. Nittler et al. 
(2001) derived "mineral mixing corrections” based on comparisons 
of theoretical spectra for mixtures of the dominant individual min- 
erals found in anhydrous meteorites (olivine, pyroxene, feldspar, 
Fe-Ni metal and iron sulfide) with those for homogeneous mix- 
tures of the same bulk compositions. The average correction factors 
(from Table 2 in Nittler et al„ 2001) are reprinted here in Table 9. If 
applied, the largest effects are to increase the S/Si and Fe/Si ratios 
by 30-40% and decrease Al/Si by 25%. We note that the magnitude 
of the Fe/Si correction factor depends strongly on the amount of 
metallic Fe present in the material. 

Recently, Foley et al. (2006) have questioned the necessity of 
the mineral-mixing correction for NEAR XRS data on the grounds 
that impact processing is likely to have reduced Eros' regolith to 
a well-mixed homogeneous powder. Moreover, calibration data for 
the Mars Pathfinder alpha-proton X-ray spectrometer (APXS) for 
powdered samples (with grain sizes of tens of microns) of me- 
teorites with a large variation in metal abundances gave correct 
compositions without requiring such a correction (Foley, 2002). 
Additional experiments are needed to clarify whether mineral- 
mixing corrections are required and we therefore report both 
corrected and uncorrected versions of the asteroid composition 
(Table 9). 

6, Fitting results for specific solar flares 

The fitting results for the eight solar flares are summarized in 
Table 8 and discussed below. Results for five of these (4-May-2000, 
15 June, 10 July, 27 December, and 2 -January-2001) are revisions of 
those previously presented by Nittler et al. (2001), For these, the 
new photon ratios remain generally consistent with published re- 
sults (Nittler et ai„ 2001, Table 3). Three additional flares for which 
data have not previously been reported (20 May, 17 June, 28 De- 
cember) are also presented here. 

Note that for many of the flare data, derived Al photon counts 
and Al/Si photon ratios were very sensitive to any limits placed on 
the zero point of the Al-fi!tered detector. For example, for the 4- 
Mav-2000 flare, calibration-source data from April 26 and May 10 
(Table 5) suggest that the Al detector zero should be above -0,04. 
When so limited, the Al/Si photon ratio is 0.092 ± 0.032. How- 
ever, if the Al-filtered detector zero is allowed to float without 
a limit, its value drops to —0,09 and the Al/Si photon ratio to 
0.066 ±0.024. The discrepancy between these solutions has been 
treated as a systematic uncertainty and is included in the total un- 
certainty quoted for this and other flares in Table 8. The effects on 
the other element ratios were minor compared to the uncertainties 
from counting statistics. 

6.1 4 May 2000 

The May 4, 2000 solar flare induced more asteroid fluorescence 
from the high-Z elements (Fe, Ca, and S) than any other solar flare 
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Meteorite Mg/Si Meteorite AS/Si 



Fig, 6, Theoretical calibration curves based on a range of bulk meteorite compositions (diamond symbols; wittier er aL 2004) and the solar spectrum and geometry appropriate 
for the solar Hare observed on 1? June 2000. Curves are quadratic fits to the data and allow the derivation of elemental abundance ratios from X-ray line flux ratios observed 
from 433 Eros. 


measured by NEAR. The measurements of the Fe/Si, Ca/Ss, and S/Si 
fluorescent photon ratios, therefore, are subject to lower statisti- 
cal uncertainties than for any other flare. However, because of the 
high off-normal angle ("panel angle"; of the gas solar monitor with 


respect to the Sun (f? meall — 35.5 s ), the derivation of the solar spec- 
trum for this flare (and therefore the conversion of the photon 
ratios from this flare to elemental abundance ratios) is also the 
most sensitive to the graded filter dimensions, 
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Table 8 

Photon counts and ratios for eight soiar flares. 
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Table 9 


Composition results: Summary (average of data from eight solar-flares; 
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The solar monitor data for this flare and a representative mode) 
(here based on Model 5, Table 3) are plotted as Fig. 7a. Because the 
Sow-energy solar flux is almost entirety cut off from the pinhole in 
the thicker-Delrin model, the emission measure of the low-T com- 
ponent must be several times higher in order to account for the 
measured low-T flux than in the thin-Delrin model. Consequently, 
the thin-Deirin filter models imply a much harder solar X-ray spec- 
trum than the thick-Delrin versions. This has a profound effect on 
the determination of the Fe/Si and Ca/Si ratios from the asteroid. 
If the thin-Delrin filter models are correct, the solar spectrum dur- 
ing this flare was quite hard, making it possible for a relatively low 
abundance of asteroidal iron to produce ail the iron signal seen in 
the asteroid-pointing detectors, in this case, the Fe/Si abundance 
implied by the asteroid fluorescence would have been relativefy 
low (Fe/Si =2.07*0.11 for Model 2; Table 11), On the other hand, 
if the thick-Deirin filter model is correct, the solar spectrum of this 
flare was much softer, requiring a higher abundance of iron in the 
asteroid in order for a relatively low number of hard solar X-rays 
to produce the iron signal detected by the NEAR XRS. Therefore, in 
this case the Fe/Si abundance ratio would have had to have been 
much higher: Fe/Si = 4,83 ± 0,27 for solar model 11. Ultimately, 
this means that the Fe/Si (and Ca/Si) results from this flare suffer 
from a large systematic uncertainty stemming entirely from uncer- 
tainty in the properties of the solar monitor graded filter, in spite 
of the good raw photon statistics from Fe and Ca. 

6.2. 20 May 2000 

The high- temperature component of the flare from May 20, 
2000 (Fig. 7b) was relatively cool and weak compared with that of 
the other studied flares. This "softer" X-ray spectrum induced few 
iron counts from the asteroid, so the statistical uncertainty in the 
derived Fe/Si value is high. A favorable spacecraft geometry with 
respect to the asteroid, however, enabled significant Mg, Al, and Si 


fluorescence to be measured by the asteroid-pointing detectors. As 
with the May 4, 2000 flare, a systematic uncertainty in the Al/Si 
photon ratio due to uncertainty in the Al-fittered detector zero is 
included in Table 8, 

6.3. 15 June and 17 June 2000 

Solar monitor data and corresponding model spectra from flares 
on June 15 and 17 of 2000 are shown in Figs. 7c and 7d, respec- 
tively. These two flares were very similar in both inferred solar 
temperature and geometry. The last calibration-source data be- 
tween June and August 2000 were taken on 8 June 2000, so the 
calibration-based detector characteristics for these flares and the 
subsequent one on July 10 are increasingly uncertain. Based on 
the unfiltered detector solutions for these flares (Table 6) it would 
appear that no resolution-degrading event affected tire asteroid- 
pointing detectors during this time period. Upper limits of 0,381 
and 0.396 were assigned to the FWHM (resolution) parameters of 
the Mg- and Ai-filtered detectors based on the 6-8 June calibration 
data. 

As with the previous flares, the systematic uncertainty in the 
Al/Si photon ratio due to Al-filtered zero calibration has been in- 
cluded in Table 8. For the June 15 flare, the high (limited AS- 
detector zero) and low (floating zero point) values for Al/Si are 
0.096 ± 0.058 and 0.054 ± 0,064, respectively. For the June 17 
flare, the high- and low-Al/Si solutions are 0.046 ± 0.060 and 
0.022*0.065, respectively. 

6.4. 10 July 2000 

Solar monitor data and a model spectrum from the solar flare 
of July 10, 2000 are plotted as Fig. 7e. Although the solar pane! 
angle for this flare was the lowest (closest to normal incidence) of 
those among the eight major solar flares, the x 2 values of the solar 
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Composition results: 
has been applied 

Effect of solar models on elemental 

ratios for 433 Eras. Each composition is the average of results from 

eight solar Hares, 

Mo mineral-mixing correction 
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Table 11 

Composition resets; individual solar flares. The model Humberts) of the graded filter model (Table 3) is given in parentheses following the corresponding ratio. 
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spectra! fitting remain relatively high regardless of the graded filter 
characteristics. The x 2 improves with degraded gas solar monitor 
resolution, suggesting that this detector’s resolution had been de- 
graded somewhat from its near-optimal values during the previous 
solar flares. Fortunately, the effect of degraded solar monitor reso- 
lution on the geochemical results is much less than the uncertainty 
due to the statistical errors in the photon counts. Since no 55 Fe cal- 
ibration data had been taken for over a month at the time of this 
flare, no limits on the zero points of the asteroid-pointing detec- 
tors were applied. The best-fit flare-derived resolution parameters 
of the asteroid-pointing detectors were 0.360, 0.373, and 0.439 for 
the unfiltered, Mg-filtered, and Aj-fiJtered detectors respectively. 
Applying an upper limit of 0.390 to the Ai-filtered detector resolu- 
tion (on the assumption that it usually did not stray too far from 


the resolution of the unfiltered detector) did not materially change 
the photon ratios. 

6.5. 27 December 2000 

This flare was unusual in that its solar monitor data could not 
be fit satisfactorily with any of the usual sets of solar elemental 
abundances (Meyer, 1985; Feldman et a!„ 1992, or the "Coronal 1" 
and ‘‘Coronal 2” sets that produced good matches to numerous 
other NEAR flares; see Table 1). The 6.5 keV iron emission structure 
was less prominent with respect to the bremsstrahlung continuum 
than was permitted by these models. Additional sets of abundances 
were tested: some with overall metallicity decreased further below 
the estimate of Meyer (1985), so that the relationships between ei- 
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ements other than H and He were unchanged; and some with the 
iron abundance reduced alone. The best fits were produced with 
only the Fe abundance reduced, and the best-performing abun- 
dance model out of the reduced-iron set (Fe/H = 1.95 x 10 -5 , as 
opposed to the Meyer (1985} value Fe/H = Si/H = 3.9 x 10"" 5 ) was 
adopted. This model is "Coronal 3" in Table 1. We note that flare- 
to- flare variability of at least this magnitude in both the Fe/Ca ratio 
and the Fe/H ratio has previously been observed (Antonucci and 
Martin, 1995; Fiudra and Schmelz, 1999). The actual and modelled 
solar monitor data are plotted in Fig. 7f, 

The solar monitor gain during this and all subsequent flares 
was 0.0415 keV/channel, an increase from the 0.040 keV/thannel 
that had consistently been observed during the summer flares. 
This may have been caused by space charge effects or other con- 
sequences of detector aging. The asteroid-pointing detectors were 
unaffected, probably because they had counted far fewer X-rays 
during the course of the mission, 

6.6. 28 December 2000 

The December 28 flare was not analyzed by Wittier et al. (2001), 
largely because the solar spectrum could not be modelled sat- 
isfactorily at that time, The CHIANTI 5.2 code yields reasonable 
models with both the “Coronal 1" and “Coronal 2" coronal abun- 
dance models. 

However, the modelling (Fig. 7g) remains less successful than 
with most of the other flares, largely because of a very large ex- 
cess of counts at about 3,5 keV. Because of the attenuation of 
the low-energy photons by the graded filter, the gas solar moni- 
tor is vastly more sensitive at 3.5 keV than it is at 2 keV. Thus, 
a large number of measured counts at 3.5 keV actually represents 
far fewer incident solar photons than a corresponding number at 
2 keV, and affects the resulting asteroid composition relatively lit- 
tle. For this reason, while fitting this flare, the region between 2.72 
and 3.44 keV was eliminated from the x 2 calculation altogether, so 
that the fits in the 2 keV and 6.5 keV regions that are more impor- 
tant to X-ray production in the asteroid would be undistorted. The 
origin of the excess signal at 3.5 keV remains unknown, though we 
note it cannot easily be explained by an increase in the solar K or 
Ar abundances in the flaring plasma, as the required enrichment 
factors (>60x for K, >5x for Ar) are implausible. 

The photon ratios in the asteroid-pointing detectors for this 
flare, unlike those for the May 4 flare, are insensitive to restric- 
tions on the asteroid-pointing detector parameters. 

6.7, 2 January 2001 

The solar spectra from the 2 January 2001 flare were easily 
fitted with the coronal-abundance solar spectral models, lacking 
either the 3 keV excess of the December 28 flare or the deficit of 
iron emission in the December 27 flare. 

The solar monitor gain during this flare was 0,0415 keV/channel. 
The FWHM parameter of the gas solar monitor was increased by 
8%. {Fitting with the degraded resolution reduced the x 2 of the 
solar monitor fitting, but the geochemical results were unchanged 
when the solar fits were done with the standard resolution (FWHM 
parameter = 0.340).) The integrated data from and models of gas 
solar monitor output are plotted as Fig. 7h. 

7. Results and discussion 

Table 10 illustrates the effect of the uncertainties in the solar 
X-ray spectra on the derived composition of 433 Eros. For each 
solar model. Table 10 lists the weighted average of the composition 
ratios derived from aii eight solar flares. These element ratios are 
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Fig. 8. XRS-derived elemental ratios for 433 Bros in the context of bulk composi- 
tions of various meteorite classes. Ellipses represent 3 o uncertainties in average 
compositions as calculated using individual solar models {Tables 3 and iO). Mo 
mineral-mixing correction (of. Foley et al., 2006; Nittler et ai, r 2G01 j has been ap- 
plied. Uncertainties in the final results {large crosshairs) combine uncertainties in 
the solar model with those due to photon statistics from the asteroid's fluorescence. 
These data are tabulated in Table 9, Top; Mg/Si vs Fe/SL Middle: Al/Si vs Ca/Si. Bot- 
tom: Mg/Si vs S/S L 


represented by the 3cr error ellipses In Fig. 8. No correction for 
mineral-mixing effects has been applied to these values. 

Results for each solar flare are individually represented in Ta- 
ble 11. For each elemental ratio derived from each flare, the maxi- 
mum and minimum values from among the solar models are tab- 
ulated. The average value for that flare from ail the solar models 
follows. The systematic errors due to solar uncertainties are es- 
timated from the maximum and minimum values derived from 





£&men;ai composi'Eio/i of 433 Eros 


143 


* Dsogetme* 

$ Howarefttes 

O Eacrites 

* H 
x L 
+ LL 
C E 
i* CV 
► CO 
$ R. 

£ AcSfRtlcttflC* 

4 Lodfaitites 
vf Winooaiies 
«■ Bradiiniit-^ 
o Ureifitss 
■q AubHics 

* Angnlcs 
$ Paibsitss 
« Mcwfeiderircs 

0.0 0J 0.2 0.3 0.4 

Al/Si 

(a) 




* Dkjgenfers 
«. Howardites 
O Eucritts 

& H 

X L 

+ U. 

D E 
o cv 
► CO 
K R 

(is Acapalcoites- 
jt, Uidranifes 
9 VVinosiaiies 
v Bradiiniffii 
c? Ureildes 
0 Aubrittvs 
a Aogritcs 
a Pslbsiws 

* M&wsidcriiss 


Fe/Si 

(b) 


Fig. 9. Flare-to-flare variation in XRS-derived elemental ratios for 433 Eros ir. the context of bulk meteorite compositions. Error bars represent to uncertainties in each of 
the eight single-flare compositions, including uncertainties due to uncertainties in the solar spectra. These data are also tabulated in Table It. 


the various solar models; these have been combined in a root- 
sumtned-square sense with the independent random errors due to 
photon statistics. (The standard deviation would not be an accu- 
rate reflection of the uncertainties due to the solar uncertainties 
because the set of models included is not a random sample of ail 
plausible solar models.) 

Table 9 presents the overall composition results of this work, 
averaged over all eight solar flares and including errors due to 
uncertainties in the solar flare spectra. The first column ("un- 
corrected") lists the compositions without any correction applied 
for mineral-mixing effects. These data are also represented as the 
large crosshairs in Fig. 8. The final column has had the mineral- 
mixing corrections applied as in Nittler et at (2001} and may 
therefore be directly compared with the results from that pa- 
per, For comparison, the average compositions of the H, L, and LL 
chondrites from the compilation of Nittler et at (2004) are also 
shown. 

Plotted in Fig, 9 are the Mg/Si, Al/Si, and Fe/Si ratios calculated 
from the eight individual solar flares, illustrating the flare-to-flare 
variability of the data. Each set of crosshairs represents the 1 <r un- 
certainty in composition for a single flare, including both errors 
due to counting statistics and those propagated from the uncer- 
tainties in the solar spectra. Although the range of element ratios 
determined from the different flares at first glance might appear to 
indicate compositional heterogeneity for Eros, in fact no one flare 


is more than 3rr away from each other or the ordinary chondrite 
field in Mg/Si or Al/Si. The flare-to-flare variability of Fe/Si is larger 
but, as discussed below, this ratio is subject to additional system- 
atic errors not related to solar modeling. Thus, the XRS data alone 
can neither establish nor eliminate the possibility of spatial varia- 
tion on 433 Eros. 

However, the NEAR MSi/NIS investigation (Veverka et al., 2000; 
Bell et ai„ 2002) has provided strong evidence for compositional 
homogeneity over ail of Eros, and the chondritic Cr/Fe ratio that 
has been derived from the XRS data (Foley et al, 2006) also sup- 
ports the hypothesis that Eros is an undifferentiated object with 
a chemical affinity to the ordinary chondrites. The averaged com- 
position in Fig, 8, therefore, probably represents the light-element 
composition of the whole of Eros. 

The composition ratios derived from individual flares have 
changed somewhat since the XRS analysis of Nittler et at. ( 2001 ), 
but the overall averages are mostly similar to the "Best Eros" 
values quoted by Nittler et al. Of particular importance, the sul- 
fur/silicon ratio remains substantially below those of chondritic 
meteorites. Troinbka et al. (2000) and Nittler et al (2001) sug- 
gested that the sulfur deficiency could be caused by either space 
weathering processes (e.g., sputtering, impact vaporization) or par- 
tial differentiation with loss of a S-rich melt. Foley et al. (2006) 
determined a chondritic Cr/Fe ratio for Eros based on XRS data 
from two solar flares and argued on the basis of this that space 
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Fig. to. Variation of Fe/Si and Ca/Si with observing geometry, sec i is the average 
secant of the incidence (off-normal] angle of the illuminarion of the portion of 
433 Eros within rhe field of view of the NtAK XR5 during each soiar flare; sec e 
is the corresponding emission angle. The data point with the highest Fe/Si ratio and 
highest (sec i -F sec e) corresponds to the solar flare of 28 December 2000: however 
the Ca/Si derived from the same flare is chondritic. 

weathering must be responsible for the sulfur deficiency since par- 
tial melting should also deplete chromium. Optical space weath- 
ering (Pieters et al., 2000; Hapke, 2001; Sasaki et ah, 2001) is 
caused by nanophase reduced iron production and deposition on 
grain surfaces. The MS1/NIS measurements (Clark et a!., 2001} es- 
tablished that the spectrum of Eros has been altered by space 
weathering. The same processes responsible for nanophase iron 
production, such as sputtering and impact vaporization, likely also 
vaporize sulfur and cause sulfur depletion in the surface lay- 
ers (e.g.. Niftier et ai„ 2001; Kilten, 2003; Kracher and Sears, 
2005). 

The new average Fe/Si ratio is slightly higher than that of 
Nittler et al. (2001), but most H chondrites fail within its Icr error 
ellipse. As discussed previously (e.g., McCoy et al„ 2001), how- 
ever, the NEAR gamma-ray spectrometer (GRS)-derived Fe/Si ratio 
of 0.8 ±0.3 (Evans et aL 2001) is considerably lower than the 
XRS derived Fe/Si and more consistent with the low-iron ordinary 
chondrite groups L and LL. Since the GRS samples to depths of 
tens of centimeters, whereas the X-rays originate from the top tens 
to hundreds of microns, this difference suggests that the elevated 
Fe/Si is a surface phenomenon. However, there are several addi- 
tional difficulties beyond solar modeling in determining Eros' Fe/Si 
ratio from XRS data, as discussed below. 

Okada (2002, 2004) has suggested that the large phase angle of 
the NEAR XRS observations should lead to erroneously high Fe/Si 
ratios due to shadowing effects. For example, in laboratory exper- 




iments, Okada (2002, 2004) demonstrated that for ground basalt 
samples sieved to various particle sizes, the ratios of Ca/Fe and 
Ti/Fe K lines decrease significantly as the phase angle of obser- 
vation increases. These authors also demonstrated that this effect 
becomes much more important as the size of the sieved particles 
increases. Okada (2002) concluded that the magnitude of this an- 
gular effect increases with the energy difference between the K 
lines that are being compared and the effect on Si/Fe (Si K a = 
1.74 keV) should thus be even more pronounced than the observed 
effects on Ca/Fe and Ti/Fe. 

Fig. 10 illustrates an attempt to explore possible geometric ef- 
fects on the Fe/Si and Ca/Si results, not through phase angle (since 
ail the NEAR data had approximately the same phase angle) but 
by plotting composition results from individual solar flares against 
a "path length parameter” which is a function of incidence and 
emission angle. This parameter, which is the sum of the secant of 
the incidence angle and the secant of the emission angle, repre- 
sents the path traveled by incoming and outgoing X-rays through 
the overburden of soil above the fluorescing atom. One would ex- 
pect matrix effects to increase as this parameter increases. The 
figure shows that the two highest single-flare Fe/Si values do come 
from the flares with the highest values of the path length pa- 
rameter. However, the Ca/Si result from the highest path-length 
flare (23 December) is apparently not inflated by this geometric 
effect 

Another effect that may be important in interpreting the NEAR 
data is size sorting in the regolith due either to the "Brazil nut” ef- 
fect (Rosato et al„ 1987) or to the reverse Brazil nut effect (Hong et 
al„ 2001). If the mineral grains in which Fe is carried have differ- 
ent physical properties than the low-Fe regolith grains (e.g,, metal 
versus silicate), size sorting in the regolith could alter the Fe/Si at 
the surface. 

Thus, given the dependence of the measured Fe/Si ratio on un- 
certainties in the solar spectrum, on phase angle effects and on the 
possibility of surface or regolith processes, the f!are-to-flare vari- 
ability in Fe/Si shown in Fig. 9 is much mote likely a reflection of 
systematic errors rather than of true elemental heterogeneity on 
the asteroid. As discussed earlier, the observed spatial homogene- 
ity in surface colors and mineralogy (based on MS1/N1S data) and 
the light-element XRS results all point to Eros being an undiffer- 
entiated asteroid. 

Overall, the two-temperature model was successful in model- 
ing the 1-10 keV NEAR solar X-ray spectra. This approach may be 
useful in the analysis of data from the Hayabusa XRS (Okada et al., 
2006), particularly in the determination of the S/Si ratio of 25143 
itokawa (or an upper limit thereof) since Hayabusa's onboard stan- 
dard sample plate lacked sulfur (Okada et al., 2006). This will en- 
able comparison of Itokawa’s S/Si value with Eros' highly depleted 
S/Si. This approach may also be helpful in analysis of Mercury data 
from the MESSENGER XRS (X-ray spectrometer; Gold et at., 2001) 
which is an improved version of the NEAR XRS design. Although 
the energy resolution of the Si-PIN solar monitor carried by the 
MESSENGER XRS is significantly improved from that of the NEAR 
gas counter, it is still insufficient to resolve the closely spaced solar 
emission lines in the region of the K-edges of the light elements 
Mg, Al, Si, and S. Thus, it will be necessary to use CHIANTI or a 
similar method of modeling the solar X-ray spectrum in conjunc- 
tion with the solar monitor data in order to obtain quantitative 
compositional results from the MESSENGER XRS, 
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